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We study a natural construction of a general class of inhomogeneous quantum walks (namely 
walks whose transition probabilities depend on position). Within the class we analyze walks that 
are periodic in position and show that, depending on the period, such walks can be bounded or 
unbounded in time; in the latter case we analyze the asymptotic speed. We compare the construction 
to others in the existing literature. As an example we give a quantum version of a non- irreducible 
classical walk; the Polya Urn. 



I. INTRODUCTION 

The study of classical random walks on a lattice has a 
long history and numerous applications in fields such as 
simulation of physical processes and probability theory. 
A random walk starts with a particle at a node on a lat- 
tice, then at each time step the particle jumps to another 
node with a given probability. Inhomogeneous random 
walks are those that that have position and direction de- 
pendent probabilities of jumping to neighbouring nodes. 

Quantum walks 0, i, S, H have been shown to have 
many applications in quantum algorithms, such as al- 
gorithms for searching for the element distinctness 
problem [6*] , for matrix product verification [7] , for test- 
ing group commutativity ^81] and for triangle finding j|. 
Quantum interference causes quantum walks to behave 
in a qualitatively different way from their classical coun- 
terparts. Most of the original work considered homoge- 
neous walks (where the amplitude for moving does not 
depend on position). In this paper we analyze a natural 
construction of inhomogeneous walks and, for a simple 
family of such walks, show that they can be bounded or 
unbounded in time. We have a number of motivations for 
this work: to gain an understanding of the most appro- 
priate general setting for quantum walks, to probe the 
possible long time behaviours of such walks, and also the 
longer term goal of producing techniques that may be 
useful for producing quantum algorithms. 

The idea of looking at inhomogeneous walks is not 
new. For example the recognition that it is natural to 
allow coins to be position dependent may be found in 
[13) [m [13] ■ alternative general method is given 
for quantizing a classical Markov chain. The method 
(particularly as generalised in [l3|) gives a large class of 
inhomogeneous walks. However as we shall see later, the 
construction in this paper is not equivalent to that in 
[isl . [1^ . As an example we give a quantum analogue 



of a reinforced process: the Polya Urn. In much of the 
earlier part of the paper we focus on walks on the line, 
but it should be clear (and the quantum Polya Urn is an 
example) that our construction is applicable to general 
graphs. 

For comparison with what comes later, we briefly re- 
view a simple homogeneous walk on the line. Let : 
n G Z} be an infinite set of states where each represents 
the position n along the infinite line of integers. In addi- 
tion there is a coin register that is spanned by two states 
(|L) and \R)) which represents the direction of motion 
at a particular time step. The evolution of the quantum 
walk proceeds by applying a coin operator to the coin 
states, in order to select which direction to move in with 
a certain probability amplitude, and then a shift oper- 
ator to move the resulting amplitudes along the line to 
their new position. The most commonly used coin op- 
erator C is the Hadamard operator, corresponding to a 
50% chance of moving left or right, and in this basis is 
given by: 



C 



\R)]{L\ 



\L) + \R)]{R\ (1) 
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The most commonly used shift operator S, which we will 
be used throughout this paper is: 

S ^^\n-l,L){n,L\ + \n+l,R){n,R\ (2) 



A full time step W of this quantum random walk is there- 
fore: 

W = S{I(S C) (3) 

So the state Y^) of the walk after T time steps is therefore: 

\^{T)) = W^\^m (4) 

A celebrated result [l^, [13] is that this homogeneous 
quantum walk spreads linearly in time, quadratically 
faster than classical walks. The shape of the distribu- 
tion is also very different, as can be seen in Figure [T] 
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FIG. 1: Overlay of distributions of classical (dotted black line) 
and quantum (solid blue line) random walks running for 200 
time steps. 

II. INHOMOGENEOUS QUANTUM WALKS 

Inliomogeneous quantum walks differ from the ones de- 
scribed above in that we allow the coin operator to de- 
pend on the position register and the coin register, rather 
than only on the coin register of the state space. This 
leads to many interesting new behaviours, such as walks 
that remain bounded in a certain region for all time, as 
we shall see. 

We will define an inhomogeneous quantum walk in a 
similar manner to the standard quantum walk, however 
we now allow the coin operator C to be dependent on m, 
the current position of the walk: 

W ^s(^\m){m\®C^ (5) 

In the case of walks on the line Cm could be an arbi- 
trary unitary operator on the two-dimensional coin space. 
Indeed, there is no need to restrict to walks in which 
one only moves to the nearest neighbours: more gener- 
ally, one could allow there to be transitions from a given 
point to any other point on the line. However for the pur- 
poses of most of our discussion of walks on the line we 
will focus on the simplest case of a transitions to nearest 
neighbours. 

As a first step in this process we will consider a family 
of examples in which the coin operator is periodic, with 
the coin operator whose matrix in the basis {\L), \R)} 
given by: 

/cos(^) -sin(if)\ 

I sin (If) cos (If) ; ^"^^ 




FIG. 2: Standard deviation of quantum walk with coin from 
equation (|6]), initial state -^(10, L) + 10, i?)) and A: = 3 against 
time steps. 

Here k is an arbitrary positive integer; the period of this 
walk is 2k. More explicitly, this means that this coin 
operator transforms states in the following way: 

— j |7i, L) + sin 1^— j |n, i?) (7) 

/ Till \ / TITT \ 

Cn\n,R) ^-sini^— j\n,L)+ COS i^—j\n,R) (8) 

Plots of the standard deviation against the number of 
time steps for the coin in equation ^ can be seen in 
figures [H [3] and [H They provide some insight into the 
behaviour that we might expect from these walks, figures 
[2] and [4] appear to be moving linearly in time and figure [3] 
appears to be periodic in time. In this section we analyze 
these cases in detail. 




FIG. 3: Standard deviation of quantum walk with coin from 
equation © , initial state ( |0, L) -|- jO, 7?) ) and k — 4 against 
time steps. 

Firstly we consider the general conditions for a walk 
to be restricted to a finite region of the line for all time: 
Lemma 2.1 Let the initial state of the walk to be of the 
form 

\no){ao\L)+bo\R)); (9) 
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FIG. 4: Standard deviation of quantum walk with coin from 



equation (|6]), initial state 
time steps. 
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(|0, L) + 10, R)) and A: = 5 against 



then an inhomogeneous quantum random walk with a uni- 
tary coin is bounded iff there exists to the left and the 
right of uq a coin with matrix of the form 







(10) 



Proof: First assume that the walk is bounded. Then 
there exists a least upper bound and a greatest lower 
bound for the walk. Let ni denote the position of the 
greatest lower bound and n„ denote the position of the 
least upper bound. Since n; is the greatest lower bound 
for the walk, the walk must reach this point and there is 
zero probability of moving further left. Hence the coin 
at ni must be of the form: 



a 
b c 



(11) 



Since the coin is unitary, it must mean that the bottom 
right component is zero as well. Thus the matrix must 
be of the following form on the lower bound: 



C„ 



e"^ 

e"t' 



(12) 



Analogously, since n„ is the least upper bound for the 
walk, it must mean that at this point there is zero prob- 
ability of moving to the right. It follows that the coin at 
n„ must be of the form: 



Therefore, if a random walk is bounded, it must have 
coins of the above form both above and below the initial 
state. 

Conversely, we now assume that there exists a coin 
with a matrix of the form (fT^ both to the left and the 
right of the initial state. Firstly consider the coin with 
form (fT2|) to the right of the initial state. The first time 
the walk approaches this position on the line, it must 
come from the left. However, the coin above makes ev- 
erything that is travelling to the right change direction 
and travel to the left. Thus if amplitude approaches from 
the left, it will never continue on to the right past that 
point. Hence this must be an upper bound for the walk. 

Similarly, if we consider the coin with the above form 
to the left of the initial state, it must first be approached 
from the right. Since the form of the above matrix flips 
the direction of movement, everything approaching from 
the right will be reflected. This means that this position 
on the line must be a lower bound for the walk. Hence, 
the quantum walk is bounded. □ 

Note that a general initial state is a superposition of 
states of the type (0) ; then the condition for boundedncss 
applies separately to each term in the superposition (so, 
in general, one could have parts of the wave- function that 
remain bounded in a region and other parts that escape). 

As an example, applying this lemma to the coin we 
deflned in equation ©, we can see that the walk will be 
bounded for even k and unbounded for odd k. 

Now that we understand when a walk will be bounded, 
we can move on to explore the behaviour when it is un- 
bounded. We would like to calculate the standard devia- 
tion of a general periodic walk in order to work out how 
fast it spreads. We are particularly interested in the fam- 
ily of walks with coin matrix © . These have even period 
for any integer k. Initially we set up the problem for gen- 
eral walks with even period 2 A. Thus we define a general 
periodic unitary coin C„.2A with period 2A, A G Z as 
follows: 



(15) 



Where C„+2A,a = Cn,A and |a„,Ap + |/3n,Ap = 1 for aU 
n. Henceforth we shall omit the A in the labelling of C, 
a and [3 in order to keep the notation tidy, i.e C„^a = Cn- 



The Period 6 Case 
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(13) 



Using the same argument as for the lower bound, the fact 
that the coin is unitary means that the coin at the upper 
bound must have the following form: 



Cn^ — 





J4>' 



(14) 



In order to illustrate the method clearly, we will first 
use the technique for the case where the period is 6, i.e. 
where A = 3. We will then move on to define the general 
case in section IIIBI 

Starting with the definition of the coin in equation (|15p 
with 2A = 6, one can derive a recurrence relation for the 
wave function: 

^{n, t+1)^ C+_,^in ~l,t) + C-+i*(n + 1, t) (16) 
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FIG. 5: Graphical representation of the computation of the 
coefficients coj in equation ((20} 



and grouping certain terms in the sum allows us to write 
the previous equation as follows: 



^a{n,t + 6) = co,o*o(»^- +co,i^'i(n- 
+ Co,2*2(« - 1, t) + co,3*o('^, t) 
+ co,4*i(ri,i) +co,5*2(ri,i) 
+ co,6*o('^ + 



(22) 



i + 6) = ci,o*i(?i -1,0+ ci,i*2(« - 1, t) 

+ ci,4*2(ri,t)+ci,5^'o(n + l,t) 
+ ci,6*i(n + l,t) 



Where: 



-^{12,1) 



(18) 



Here tphin, t) is the amplitude of the wave function trav- 
elling left at position n and time i, and il)R{n,t) is the 
corresponding amplitude travelling right, i.e. 

oc 

m))^ {'4'L{n,t)\n,L)+xl^R{n,t)\n,R)) (19) 

Recursively substituting the recurrence equation into 
itself six times, starting from ^'(6n, 6), 5'(6n + 2,t + 6) 
and \l/(6n + 4, t + 6) yields equations of the form: 

6 

«'(6n, i + 6) = ^ co,,^'(6(n - 1) + 2j, t) 

3=0 
6 

*(6n + 2,t + 6) = Y^ cij^{6{n - 1) + 2j + 2, t) (20) 

6 

*(6n + 4, i + 6) = ^ C2 j'I'(6(n - 1) + 2j + 4, t) 

where the Cij are constant (i.e. independent of n) 2 x 
2 matrices given by sums of products of C+ and Cg . 
There is a simple way of computing these coefficients by 
referring to Figure [5] . For example 

/-^+ /-^+ /-^+ 

Co,0 — ^5 '^4 ^3 '-'2 '-'1 "-"O 

„ 1^+ /^+ r^— 

Co,l — L'4 ^3 "-^2 "-^l '-'2 

+05 O4 O3 O2 O3 O2 

-I-O5 O4 O3 O4 O3 C-2 

I /^+ /^+ 

-I-O5 O4 O5 O4 O3 C-2 

_L/^+ 

-I-O5 O5 O4 O3 C-2 

+CrC+C+C+C+C+ (21) 
We will change notation at this point to make the equa- 
tions more concise. We define 5'(6n -I- 2i,t) := '^i{n,t) 



«'2(n, t + C2,o4'2(ri - l,t) + 024*0(^1, 
+ C2,2*l("-,0 +C2,3*2(n-,i) 
+ C2,4*o(" + 1,0+ C2,5*l(»^ + 1, t) (24) 
+ C2,6*2("-+ l,i) 

Discrete Fourier transforming these equations yields the 
following: 

*o('^, i + 6) = (co,oe-''" + co,3 + co^ee*") *o('^, t) 

+ (co4e"'"+co,4)*i(w,0 (25) 
+ (00,26"*" + co^s) *2(t^,i) 

^i(cj, t -f 6) = (ci,2 + ci^se*") *o(w, 

+ (ci,oe-''" + ci,3 + ci,6e*") *i(^, t) (26) 
+ (ci,ie-*" + ci,4) *2(w,i) 

^2(W, t + 6) = (C2,l + C2,4e*") *o('^, t) 

+ (c2,2+C2,5e^")^'i(co,t) (27) 

+ (C2,0e"*" + C2,3 + 02,66*'") ^2i^, t) , 



where 



n — — 00 

1 



/(n) ^ — I dio e"^'^f{uj). (28) 



Thus we can give an expression for 4'(a;,6T -I- 7) in 
Fourier space where {7 = 0, 1, 2, 3, 4, 5}: 

^'o(w,6r + 7)\ /«'o(c^,7)' 
^'i(a;,6r-K7) = «'i(^,7) | (29) 

^'2(u;,6r + 7)/ V^'2(w,7) 
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Where the 6x6 matrix U is: 



co.oe + Co, 3 + co,6e'" co.ie 
U = I ci,2 + ci,5e^- Ci,oe-*" + 

C2,l + C2,4e'" C2,2 4 



The matrix [/ in equation (j29p , is unitary, which allows 
us to write it as: 



+ Co,4 

ci,3 + cifie" 



C0,2e + Co.5 

ci.ie-*" + ci,4 



(30) 



f/(cc;) = Uo{iu)D{cu)U^{cu) 



(31) 



where C/q is a unitary matrix composed of the eigenvec- 
tors of U and Z? is a diagonal matrix containing the eigen- 
values of U. This allows us to rewrite U'^ in the following 
way: 



(32) 



Since the eigenvalues of a unitary matrix lie on the unit 
circle on the complex plane, 



(33) 



where Sjk is the Kronecker delta. 

This yields integral expressions for the wave function 
at arbitrary time and position: 



V'ol(", 6T) 



1 



6 



27r / 

1=1 



ai,z(c^)e^^(^'(")+'^«d^, (34) 



■'-^ 1=1 



ipiL{n 



:6T) = — / VQ;3,i(w)e 



inMu.)+u0)^^^ (36) 



(„,6T) = ^a4,i(^)e^^(^'(-)+-«da;, (37) 



V'2i?(n,6T) = r V «6,z(w)e'^(^'(")+"«da;, (39) 



where: 
ay (w) : 
(f/o)y- 

and 



^(C/d)j,2m+lV';L(t^,0) + (f/d)j,2m+2V';fl,(t^,0) 

(40) 



.m=0 



/3 = n/T. 



(41) 



We now use the integral expressions in equations 
to ((39)) to calculate the first and second moments of the 
probability distribution on the line. From these we can 
calculate the standard deviation and thus the rate at 
which the walk spreads. It is important to note here 
that we compute these moments only for a number of 
time steps that is a multiple of six. 

The probability p that the walk is at position 6?! -I- 2j 
after 6T time steps is given by: 

p(6n -f 2j, QT) = |7/.,x(n, 6T)|2 + |V7,fl(«, 6T)|2 (42) 

Hence we can calculate the first moment after 6T time 
steps: 



E(iV,6r)= E(6« + 2j)p(6« + 2j,6T). (43) 

n— — oo j—Q 



,6T)^^ f y as i(w)e^'^(^'(")+"'5'dw, (38) Or more explicity, 
■^^^ 1=1 



oo 6 6 



E(iv,6r)= E E E 



n= — oo j = l li^l2 = l 



6n + 2L^J 
(27r)2 



duji / dw2e^^"^""^^"aj,ii(wi)a^Ju;2)e*^(^'i('"i)-^'2('"^» (44) 



Where is the truncated integer part of 



We are interested in the leading order of equation ([44]) 
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in T as T ^ oo because we are calculating the long time 
variance of the system. We can simplify this expression 
by noting that the integrals with coefhcient (jfp' I'^^l 
in equation (|44|) are of order 0(1) using the method of 
stationary phase. 

As we will show that 0(1) is not the leading order in 
T, we can discard it from our calculation. We can further 
simplify the expression using: 



1 



2tt ^ 

X— — 00 



oo 



')x _ 



(-z)"(5(")(a;-cj') 



(45) 



This yields a much simpler equation for the dominant 
terms in the expectation of N: 



E{N, 6T) 



3z f'^ 



(46) 

For li = I2, the exponentials disappear and we are 
left with terms linear in T or constant after the integral 
over uji has been performed. Hence these terms give a 
leading order term E,{N, 6T) ~ doT assuming that do 
does not vanish. When h ^ I2 we can apply the method 
of stationary phase as T ^ 00, to each term, assuming 
that the original U is non-degenerate (as we will find 
in the example of particular interest to us - see below) 
and find that each term is, at most, 0{^/T). Hence we 
can say that the expectation of N has leading order 



E{N, 6T) 



6^ /'^ ? E M^)\'>^i{^)- 

■^-^ 3,1 = 1 



(47) 



as T — » 00. 

The second moment for an even number of time steps 
can be calculated in a similar way using: 

2 

¥.{N\QT) = E(6'^ + 2jf p(6n + 2j, 6T) (48) 

n j=0 

However, the 24nj and 4j^ coefficients of p are at most 
0{T) and 0(1) respectively. Since we will show that 
neither of these are leading order in T we can ignore them 
from the calculation and end up with an expression for 
the leading order behaviour for the second moment: 

E(iV^6r) 

(49) 

Using the same reasoning as before, when li = I2 we 
have terms like d2T^ + d^T + ^4 where (^2,^3 and ^4 are 



constants. Using the method of stationary phase for the 
other terms yields terms of order less than or equal to 
T 2 . Hence, assuming ^2 does not vanish and the leading 
order term in E{N'^,6T) behaves like as T ^ c». 

Collecting the leading order terms together we find 
that the leading order term in the variance is 



Var(7V, 



where 



GT) ^ ^ / dcjY \adu^)\'{Ki^) ~ , 

(50) 

1^^^ r du;Y\M^)\'>^i(^)- (51) 



The expression (I50p is valid for any coin of period six. 
And given its form - an integral of a sum of terms, each 
of which is a product of two positive terms - one can 
see why one might expect, typically, that the leading or- 
der behaviour of the standard deviation will indeed to be 
linear in T. We have not been able to get closed form 
analytic expressions for aji{uj) or AJ(aj) (they involve di- 
agonalising a 6 x 6 matrix). However in the case that 
the coin Cn has the matrix with the simple form (with 
period 6) 



3 J 
n-TT \ 



COS 



A 3, 



and where the initial state is 



\m) = io) 



1^) + 1^) 

V2 



(52) 



(53) 



we give figures below plotting X)j=i |Q'j;('^)P and AJ(w) 
for / = 1...6. These show that the variance does in- 
deed have leading order behaviour proportional to T^, 
and hence the spread of the walk is proportional to T, as 
observed in the earlier figure [H 




FIG. 6: Plot of and Ylj^^ \aj2{uj)r. Graphs 

of both are identical, hence the single graph representing ei- 
ther sum. The initial state is (|0, L) + |0, R)) 



FIG. 7: Plot of Ej=i lai3(t^)|^ and ^ {ajiHf- Graphs 
of both are identical, hence the single graph representing ei- 
ther sum. The initial state is -7= (|0, L) + |0, R)) 



\ 20 - 




\ ~ 




\ 10 - 








1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



-3-2-10 1 2 



FIG. 8: Plot of Y:%i \aj5H\'^ and Yl%i Graphs 
of both are identical, hence the single graph representing ei- 
ther sum. The initial state is -j= (|0, L) -I- \0,R)) 



FIG. 9: Graphs of AJ(lj). Each line represents the graph of a 
different value of I. The black line is X',, the blue line is A2, 
the brown line is A3, the red line is A4, the green line is A5 
and the magenta line is Ag. 

the 2A X 2A unitary matrix defining the time evolution: 

/ ^'o(c^,2Ar + 7) \ / ^-0(^^,7) \ 

: Uf/H (56) 

V*A-i(w,2AT + 7) / V*A-i(w,7)/ 

Earlier it was shown that walks with coins of the form 
given in equation ([6]) with k even are bounded for all time 
(see also Fig. [3|). A natural question that arises is how 
this fits in to the derivation above. 

Obviously the "leading term" in ((54|) must be zero. 
Clearly there are various ways that this could occur; in 
particular ([54]) will be zero if the A/(w) are in fact inde- 
pendent of Lv. By direct calculation we have checked that 
this is indeed the case for k — 2 and k — A. 



B. The General Case 

Very similar calculations can be performed for the pe- 
riod 2A case, for general A, leading to an expression for 
the leading order behaviour exactly analogous to that for 
the period six case: 

Var(iV, 2AT) ~ 

3,1=1 

where 

2A 

^, = ^ r rf'^ E k-^hi'a;(^)- (55) 



III. QUANTUM WALKS DEFINED VIA 
DOUBLE- REFLECTION 

In [1^, Szegedy proposed a method for defining the 
quantum walk on a graph starting with a classical 
Markov chain. The work was generalized in [ij]- For 
reasons that will be clear shortly, we refer to these quan- 
tum walks as "double- reflection" walks. 

We now compare the construction we have been using 
with that in (isl . [l^ for the particular case of motion 
on the line. The state space used in [l^, consists of 
two copies of the position space, rather than a position 
and coin register. Thus the Hilbert space is spanned by 
|n,m), n,m = — c»...c». We can describe the coin based 
walks using this Hilbert space by identifying \n, L) with 
|n, n + 1) and |n, R) with |n, n — 1). 

In the case of motion on a line, a double-reflection 
walk is set up as follows [l^ 13 ■ Let c?„ , e„ , /„ and g„ 
be complex numbers satisfying 



The functions ajiioj) and \i{lo) arise, as in the period 
six case, from the initial condition and diagonalization of 



Mn|' + |e„P = l,V7ieZ 

|/n|' + l5«P = l,VneZ 



(57) 
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The two-step walk operator Wdr is defined via the 
foUowing equations: 



\Pn) = dn\n + 1) + e„|n - 1) 
km) = fm\m + 1) + gm\m - 1) 

ns = ^ \q7n)\m){qm\{m\ 

711 

Wdb. ^ {2Ub - I){2Ua - I) 



(58) 



(59) 



(60) 



The Hadamard walk in the double-reflection frame- 
work can be realized by setting: 



(61) 



The inhomogeneous quantum walk with coin in equation 
([6|) can be reproduced by taking: 



1 



d„ = /„^-^^l + sin(^) 
1 I 7rvn\ 



(62) 



Although the original walks in [l3| have real ampli- 
tudes, it is natural to allow the constants c?n,en,/n and 

be complex numbers as in 14]. We now show that even 
if we allow this, not all walks produced by the position- 
dependant coin construction can be realized by walks of 
the form ((60)) . 

Using a generalized double-reflection form of the quan- 
tum walk to IV'o) = OLo\n, n -I- 1) -I- 5o|n, n — 1) yields: 

W^Dfl.lV'o) = 

Vn-Wn-i {2aodnen + ^0 (2|e„p - l)) |n - 2, n - 1) 
+ (2|/„-iP - 1) (2ao4e„ + 60 (2|e„|2 _ 1)) |n, n - 1) 
+ (2|g„+i|2 - 1) (ao (2|d„|2 „ 1) ^. 26od„e„) \n,n+ 1) 

-h2/„+i3„+i (ao (2|(i„|2 - 1) + 26od„e„) |n + 2,n + 1) 

(63) 

We now wish to compare this to a walk using a con- 
trolled unitary coin. Since we have allowed the "forward" 
and "backward" steps to be different in we need to 
take this into account. So for two steps of the position- 
dependent-coin walk we use different coin operators in 
the first and second step. 

We consider a two general controlled unitary coins, de- 
scribed in the Hilbert space of two copies of the position 
register: 

Un\n,n+ 1) = an\n,n + 1) + (3n\n,n- 1) 

Un\n,n- 1) = -e'^"/3„|n,n 1) + e'^"a„|n,n - 1) 

(64) 



and 

Un\n, n -I- 1) = a„\n, n + I) + f3„\n, n - 1) 

IJn\n,n- 1) = -e'^"/3„|n,n + 1) + e*^"a„|n,n - 1) 

(65) 

The general two-step walk operator associated with 
these coins is given by: 



where 



WpDC = SCSC 



C= |m)(m|(g)J7„ 



m— — 00 



C= ^ \m){m\(g)ij„ 



(66) 



(67) 



(68) 



m— — 00 



and 



5"= |m — 1, m) (m, TO + 1| + |m + 1, TO)(m, m — 1| 



\mn){mn\ 

m^n-\-l and m^n — 1 



(69) 



Let \ipa) = ao|n, 7i+l)+6o|n, n-1) with |aop-l-|6oP = 
1. Applying walk operator yields: 

WpDc\lpo) = 

a„-i (aoa„ - boe''^"f3„) |n - 2, n - 1) 
-|-^„_i (aoa„ -5oe"'"/3„) |n,n- 1) (70) 
- e*^"+i/3„+i (ao/3„ ?>oe'''"a„) |n, n-f 1) 
+ e'"'"+i(5„+i (ao/3„ + 5oe^^"a„) |n + 2, ?i + 1) 

We can now see that the effect of Wdr. can be re- 
produced by a walk using a position-dependent coins by 
choosing a„ = 2J„e„, /3„ = - (2|e„p - l), 6'„ = and 
a„ = 2/„g„, /3„ = - (2|g„|2 - l), 0„ = 0. 

However the most general walk of the form cannot 
reproduce all walks with position dependent coins. This 
is because from (l63l) one can see that 



{n,n+l\WDR\n,n+l) = (2|5„+i|2-l) {2\dnf - l) 

(71) 

is real for any choice of d„ , e„ , /„ and gn , however 

{n,n+l\WpDc\n,n + l) = -e'^-+'Pn+iPn, (72) 

which need not be real. Indeed even if we take both steps 
of the position-dependent-coin walk to be the same (as 
is usually done), we get 

{n,n+l\WpDc\n,n+l) ^ -e'^'-^^Pn+iPn, (73) 



9 



and still the most general walk of the form ((63)) cannot 
reproduce all walks with position dependent coins. 

Indeed it is not too difhcult to check that even if one 
allows transitions to arbitrary positions in the double- 
reflection form of the walk, it cannot reproduce all walks 
defined by position-dependent coins. For let us take gen- 
eral states of the form 



\Pn) = DnM + j) 

oo 

IQm) = X! ^rnj\m+j) 
j=-oc 

rather than (1581). The matrix element 



(74) 



{n,n + l\WDB.\n,n + 1) 

= (2|(n|Q„+i)|2-l)(2|(7i + l|P„)p-l) 



(75) 



is still real. 

We also note that for some purposes it is natural to 
use the shift operator 



5 = ^ |n - 1, R}{n, L\ + \n + 1, L) (n, R\ 



(76) 



rather than ^ . This is a closer analogue of the shift em- 
bodied in the double-reflection framework. It is straight- 
forward to check, however, that this does not alter 
our conclusions: this modified version of the position- 
dependent-coin walk can also reproduce the double- 
reflection walk; and not all position-dependent-coin walks 
with this modified shift can be realized by double- 
reflection walks. 

Nonetheless, there may be a generalization of the 
double-reflection form of the quantum walk which allows 
for these two constructions to be equivalent; this question 
is left open for future work. 



IV. THE POLYA URN 

The classical Polya Urn is a model in which there is an 
urn containing a number of red and black balls. At every 
time step one is drawn at random, then replaced with 
a copy of the drawn ball. Let i?„ be a random variable 
representing the number of red balls in an urn containing 
n balls in total. Let the random variable Xn be given by 
= then it is possible to show that /3(?'o: ^o) 

almost surely as n — > cx3, where vq and bo are the initial 
number of red and black balls in the urn respectively, and 
/3(rni bo) is a beta distribution with parameters rp and bo- 

The classical system is not irreducible so it does not 
fit naturally into the double-reflection framework. It is 
nonetheless instructive to see what happens if we try to 
use that framework to quantize this system. 

A natural position space is spanned by |r, b) where r 
and b are the numbers of red and black balls respectively. 



In the classical model a time step increases the number of 
red or black balls by one so it is natural to try to set up 
a quantum model by taking something like the following: 

\Pr,b) = ar,b\r + !,&>+ Pr.b\r, b+1), 

\qr.b) = 7r,6|»' + !,&)+ Srfi\r,b+ 1), 



(77) 



Ua = Y^ \r,b){r,b\ (g) \pr,b){Pr,b\^ 

r,b 

^Yl \lr,b){qr,b\ (E) \r,b){r,b\, 

r,b 



Wf 



(2nB-/)(2n^-/), 



(78) 



(79) 



'Polya-DR 

where cnr,b^ lir,b^lr,b^5r,b are constants. 

However it may easily be checked that this system has 
the following behaviour 

• for each pair (r, b) the two-dimensional subspace 
spanned by jr, 6)|r + 1,6) and \r^b)\r,b + 1) is in- 
variant under the time evolution (so if the starting 
state of our walk is in this subspace, it just moves 
around in the subspace) 

• similarly for each pair (r, 6) the two-dimensional 
subspace spanned by |r-|-l, 6) |r, 6) and |r, 5-1-1) |r, 6) 
is invariant under the time evolution 

• all other states are unchanged by the time evolution 

In particular the number of red or black balls does 
not increase with time by more than one from its initial 
value. Thus this quantum evolution does not seem to 
be the natural quantum version of the classical model. It 
would be interesting to know whether it is possible to use 
different states than ((77|) to produce a double-reflection 
quantization of the Polya Urn that is more satisfactory 
(i.e. one that only increases the number of red and black 
balls as time evolves), but we have not yet been able to 
do so. 

It turns out to be relatively straightforward to set up 
a quantum walk with a position-dependent coin that 
seems to capture how a quantum version of the Polya 
Urn should behave, as we now describe. 

The state space is spanned by |r, 6, i?) and |r, 6, B). 
The first two registers are the number of red and black 
balls and the third register can be in one of two states 
R or B corresponding to whether the number of red or 
black balls is to increase in the "shift" step of the walk. 
We are really only interested in the number of red and 
black balls r and b being positive, but in order for the 
walk to be unitary it will be convenient to allow r and b 
to take any integer values (although we will typically only 
be interested in initial states with r, 6 > and r-\-b > 1). 
So we set up our walk as follows. The coin operator is 
defined by: 



Cr^b\R) 



'7T-b\R) + j7T-b\B) ifr,6,r + 6-l>0. 



otherwise. 



(80) 
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\\B) 



othcrwist 



5*= J2 \r+l,b,R){r,b,R\ + \r,b+l,l 



r,b—~oo 



Wpolya-PDC = 

\ r,6— — oo 

An attractive feature of this walk is that 
one step of the walk, measure the system and reset the 
chirality state to \R), then the classical walk is recovered. 

The resulting probability distribution of running this 
walk, starting with 10 red balls and 10 black balls for 200 
time steps, can be seen in Figure [TOl 




20 40 60 80 100 120 140 160 180 



FIG. 10: Probability distributions for number of red balls of 
the quantum Polya Urn (dashed red line), number of black 
balls of the quantum Polya urn (solid black line) and the 
number of red balls of the classical Polya Urn (dot-dashed 
blue line) with ro = feo = 10 for 200 time steps. 

There is a heavy bias to the right, similar to the 
Hadamard random walk. It is possible to reduce the bias 
using a more symmetric coin and initial condition. 

The quantum Polya urn plot shows several features 
that one might expect, such as most of the probability 
being concentrated to the far right or left. This is to be 
expected, as the Polya urn is an example of a reinforced 
process; if an event occurs then it becomes more likely to 
occur in the future. 

We also calculated the standard deviation of the Polya 
Urn quantum walk numerically, and a plot of this divided 
by the number of timesteps can be seen in Figure 1111 
For the walk to evolve linearly in time asymptotically as 
T — > cxD, one would expect this graph to tend towards 
a horizontal line. This appears to be the case, but a 
proof of asymptotic linear evolution is an open question 
for future work as the quantum Polya Urn has a non- 
periodic coin and thus our original method for periodic 
coins can not be applied. 
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FIG. 11: Standard deviation divided by number of time steps 
for the quantum Polya Urn with vq = 10 and bo — 10 run for 
20000 time steps 



V. CONCLUSIONS 

In this paper we have analyzed some particular models 
of inhomogeneous walks. It would be interesting to un- 
derstand whether one can make general characterizations 
of the long-time behaviour of quantum walks from knowl- 
edge of their (position-dependent) coins. For example we 
have described walks that are bounded and others whose 
spread is linear in time. We do not know what interme- 
diate types of behaviour are possible. It would also be 
interesting to know whether quantum walks (including 
non-periodic ones) "typically" spread linearly in time. It 
would also be attractive to understand the relationship 
between the double-reflection and position-dependent- 
coin walks more fully: is there a suitable generalization 
that encompasses both? 
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